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Abstract 

The methods developed in the previous Part I are applied to a few impor- 
tant reductions of BBGKY hierarchy, namely, various examples of Vlasov-like 
systems. It is well known that they are important both for fusion modeling and 
for particular physical problems related to plasma/beam physics. As in part I 
we concentrate mostly on phenomena of localization and pattern formation. 
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LOCALIZATION AND FUSION MODELING 
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REDUCTIONS* 



The methods developed in the previous Part I are applied to a few 
important reductions of BBGKY hierarchy, namely, various examples of 
Vlasov-like systems. It is well known that they are important both for fu- 
sion modeling and for particular physical problems related to plasma/beam 
physics. As in part I we concentrate mostly on phenomena of localization 
and pattern formation. 



1 INTRODUCTION: VLASOV-POISSON 
SYSTEM 



1.1 Description 

In this part we present the applications of our approach based on variational 
multiresolution technique [l]-[6], considered in Part I [7], to the systems with 
collective type behaviour described by some forms of Vlasov-Poisson/Maxwell 
equations, some important reduction of general BBGKY hierarchy [8]. Such 
approach may be useful in all models in which it is possible and reasonable to 
reduce all complicated problems related to statistical distributions to the prob- 
lems described by the systems of nonlinear ordinary/partial differential/integral 
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Abstract 



equations with or without some (functional) constraints. In periodic accelera- 
tors and transport systems at the high beam currents and charge densities the 
effects of the intense self-fields, which are produced by the beam space charge 
and currents, determinine (possible) equilibrium states, stability and transport 
properties according to underlying nonlinear dynamics. The dynamics of such 
space-charge dominated high brightness beam systems can provide the under- 
standing of the instability phenomena such as emittance growth, mismatch, halo 
formation related to the complicated behaviour of underlying hidden nonlinear 
modes outside of perturbative tori- like KAM regions [8]. Our analysis based 
on the variational-wavclct approach allows to consider polynomial and rational 
type of nonlincaritics. In some sense in this particular case this approach is 
direct generalization of traditional nonlinear SF approach [8] in which weighted 
Klimontovich representation 

S fj = a 3 W H S ( X - X ji) S (P - Pji) (!) 

or self-similar decompostion like 

Srij = b j Y,w ji s(x-x ji ), (2) 

where s(x — Xji) is a shape function of distributing particles on the grids in 
configuration space, are replaced by powerful technique from local nonlinear 
harmonic analysis, based on underlying symmetries of functional space such as 
affine or more general. The solution has the multiscalc/multircsolution decom- 
position via nonlinear high-localized eigenmodes, which corresponds to the full 
multircsolution expansion in all underlying time/phase space scales. Starting 
from Vlasov-Poisson equations, we consider the approach based on multiscale 
variational-wavelet formulation. We give the explicit representation for all dy- 
namical variables in the base of compactly supported wavelets or nonlinear 
eigenmodes. Our solutions are parametrized by solutions of a number of re- 
duced algebraical problems, one from which is nonlinear with the same degree of 
nonlinearity as initial problem and the others are the linear problems which cor- 
respond to the particular method of calculations inside concrete wavelet scheme. 
Because our approach started from variational formulation we can control evo- 
lution of instability on the pure algebraical level of reduced algebraical system 
of equations. This helps to control stability/unstability scenario of evolution in 
parameter space on pure algebraical level. In all these models numerical mod- 
eling demonstrates the appearance of coherent high-localized structures and as 
a result the stable patterns formation or unstable chaotic behaviour. Analysis 
based on the non-linear Vlasov equations leads to more clear understanding of 
collective effects and nonlinear beam dynamics of high intensity beam propaga- 
tion in periodic-focusing and uniform-focusing transport systems. We consider 
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the following form of equations 
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The corresponding Hamiltonian for transverse single-particle motion is given by 



H(x,y,p x ,p y ,s) = -(p 2 x +p 2 y ) + ^[k x {s)x 2 
+k y (s)y 2 } + H 1 (x,y,p x ,p y ,s) + tp(x,y,s), 



(4) 



where Hi is nonlinear (polynomial/rational) part of the full Hamiltonian and 
corresponding characteristic equations are: 



d 2 x , . . d , , . 
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dfi2 +k y (s)y+ —tp(x,y,s) = 0. 



1.2 Multiscale Representation 

We obtain our multiscale/multiresolution representations for solutions of these 
equations via variational-wavelet approach. We decompose the solutions as 



where set 



fb(s,x,y,p x ,p v ) = ^2 ®3\f( s , x ,yiPx,P y ), 

i—i c 

oo 

ip(s,x,y) = ^2 ®5 J tp(s,x,y), 

j=jc 

OO 

x(s) = ^2 ®$ k x(s), 

oo 

y(s) = J2 ©^V*). 
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i^ci jci k Cl £ c ) 
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(7) 

corresponds to the coarsest level of resolution c in the full multiresolution de- 
composition [9] 

V c C V c+1 C V c+2 C . . . (8) 
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Introducing detail space Wj as the orthonormal complement of Vj with respect 
to 

V j+ i:V j+1 =V j @W j , (9) 

we have for 

f^,x,ycL 2 (K) (10) 



L 2 (R) = V c ($W j . (11) 



3=c 



In some sense it is some generalization of the old SF approach [8] . Let L be 
an arbitrary (non) linear differential/integral operator with matrix dimension 
d, which acts on some set of functions 

* = *(«, a;) = (* 1 (s,a;),...,* d (s,a;)), s,a;eftcR n+1 (12) 

from L 2 (Q): 

LV = L(R(s,x),s,x)V(s,x)=0, (13) 

where x are the generalized space coordinates or phase space coordinates, and 
s is "time" coordinate. After some anzatzes the main reduced problem may be 
formulated as the system of ordinary differential equations 



Qi(f)^ = Pi(f,s), / = (/l,- -,/„), (14) 
i = 1, . . . ,n, max (leg Pi = p, maxefeg Qi = q 

i ' i 

or a set of such systems corresponding to each independent coordinate in phase 
space. They have the fixed initial (or boundary) conditions fi(0), where Pi,Q% 
are not more than polynomial functions of dynamical variables fj and have 
arbitrary dependence on time. As result we have the following reduced algebraic 
system of equations on the set of unknown coefficients Af of localized eigenmode 
expansion: 

L(Q ij ,X,a I )=M(P ij ,X,p J ), (15) 

where operators L and M are algebraization of RHS and LHS of initial problem 
and A are unknowns of reduced system of algebraical equations (RSAE). After 
solution of RSAE (15) we determine the coefficients of wavelet expansion and 
therefore obtain the solution of our initial problem. It should be noted that if 
we consider only truncated expansion with N terms then we have the system of 
N x n algebraic equations with degree 

£ = max{p, q} (16) 
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and the degree of this algebraic system coincides with degree of initial differential 
system. So, we have the solution of the initial nonlinear (rational) problem in 
the form 

N 

/*(*) = . MO) (17) 
k=i 

where coefficients Xf are the roots of the corresponding reduced algebraic (poly- 
nomial) problem RSAE. Consequently, we have a parametrization of solution 
of initial problem by the solution of reduced algebraic problem. The obtained 
solutions arc given in this form, where fk {t) are basis functions obtained via mul- 
tiresolution expansions and represented by some compactly supported wavelets. 
As a result the solution of equations has the following multiscale/multiresolution 
decomposition via nonlinear high-localized eigenmodes, which corresponds to 
the full multiresolution expansion in all underlying scales starting from coarsest 
one. For 

x = {x,y,p x ,p y ) (18) 

we will have 

*(«,x) = flijU^^'^x), (19) 

y 3{s) = vi/ slow (s) + J2v i j (^), ^~2< 

1>N 

IP(x) = U*(x)+ £ l4(fc ro x), fc m ~2"\ 

m>M 

These formulae give us expansion into the slow part and fast oscillat- 

ing parts for arbitrary N, M. So, we may move from coarse scales of resolution 
to the finest one for obtaining more detailed information about our dynam- 
ical process. The first terms in the RHS correspond on the global level of 
function space decomposition to resolution space and the second ones to de- 
tail space. It should be noted that such representations give the best possible 
localization properties in the corresponding (phasc)spacc/time coordinates. In 
contrast with other approaches this formulae do not use perturbation tech- 
nique or linearization procedures. So, by using wavelet bases with their good 
(phase) space/time localization properties we can describe high- localized (co- 
herent) structures in spatially-extended stochastic systems with collective be- 
haviour. Modeling demonstrates the appearance of stable patterns formation 
from high- localized coherent structures or chaotic behaviour. On Fig. 1 we 
present contribution to the full expansion from coarsest level of decomposition. 
Fig. 2 shows the representations for full solutions, constructed from the first six 
scales (dilations) and demonstrates (meta) stable localized pattern formation 
in comparison with chaotic-like behaviour (Fig. 10, Part I) outside of KAM 
region. We can control the type of behaviour on the level of reduced algebraic 
system (15). 
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2 RATE/RMS MODELS 



2.1 Description 

In this part we consider the applications of our technique based on the meth- 
ods of local nonlinear harmonic analysis to nonlinear rms/rate equations for 
averaged quantities related to some particular case of nonlinear Vlasov-Maxwell 
equations. Our starting point is a model and approach proposed by R. C. 
Davidson e.a. [8]. We consider electrostatic approximation for a thin beam. 
This approximation is a particular important case of the general reduction from 
statistical collective description based on Vlasov- Maxwell equations to a fi- 
nite number of ordinary differential equations for the second moments related 
quantities (beam radius and emittance). In our case these reduced rms/rate 
equations also contain some distribution averaged quantities besides the second 
moments, e.g. self-field energy of the beam particles. Such model is very ef- 
ficient for analysis of many problems related to periodic focusing accelerators, 
e.g. heavy ion fusion and tritium production. So, we are interested in the under- 
standing of collective properties, nonlinear dynamics and transport processes of 
intense non-neutral beams propagating through a periodic focusing field. Our 
approach allows to consider rational type of nonlinearities in rms/rate dynami- 
cal equations containing statistically averaged quantities also. The solution has 
the multiscale/multiresolution decomposition via nonlinear high- localized cigen- 
modes (waveletons), which corresponds to the full multiresolution expansion in 
all underlying internal hidden scales. We may move from coarse scales of resolu- 
tion to the finest one to obtain more detailed information about our dynamical 
process. In this way we give contribution to our full solution from each scale of 
resolution or each time/space scale or from each nonlinear eigenmode. Starting 
from some electrostatic approximation of Vlasov-Maxwell system and rms/rate 
dynamical models we consider the approach based on variational-wavelet for- 
mulation. We give explicit representation for all dynamical variables in the 
bases of compactly supported wavelets or nonlinear eigenmodes. Our solutions 
are parametrized by the solutions of a number of reduced standard algebraic 
problems. 

2.2 Rate Equations 

In thin-beam approximation with negligibly small spread in axial momentum 
for beam particles we have in Larmor frame the following electrostatic approx- 
imation for Vlasov-Maxwell equations [8] : 






(21) 
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where ip{x, y, s) is normalized electrostatic potential and F{x, y, x' , y' ', s) is dis- 
tribution function in transverse phase space (x, y, x', y', s) with normalization 

N = J dxdyn, n(x,y,s)= Jdx'dy'F, (22) 

where K is self-field perveance which measures self-field intensity. Introducing 
self-field energy 

E ( s ) = / dxdy\d 2 ?P/dx 2 +d 2 ?P/dy 2 \ (23) 

4ttK J 

we have obvious equations for root-mean-square beam radius R(s) 

R(s)=<x 2 +y 2 > x ' 2 (24) 

and unnormalized beam emittance 

e 2 (s) = 4(< x' 2 + y' 2 >< x 2 + y 2 > - < xx' - yy' >), (25) 

which appear after averaging second-moments quantities regarding distribution 
function F [8]: 

d 2 R(s) /,,.„, . AYl + A)\„, , e 2 (s) 

dg2(s) i °R 2 (,\( dR K{1 + A) dE{s) \ n 
__ + 8i?(s) (___ _j _ o, 

where the term 

Jf(l + A)/2 (27) 
may be fixed in some interesting cases, but generally we have it only as average 

K(l + A)/2 = - < xdtp/dx + yd7p/dy> (28) 

w.r.t. distribution F. Anyway, the rate equations represent reasonable reduc- 
tions for the second-moments related quantities from the full nonlinear Vlasov- 
Poisson system. For trivial distributions Davidson e.a. [8] found additional 
reductions. For KV distribution (step- function density) the second rate equa- 
tion is trivial, 

e(s) = const (29) 

and we have only one nontrivial rate equation for rms beam radius. The fixed- 
shape density profile ansatz for axisymmetric distributions also leads to similar 
situation: emittance conservation and the same envelope equation with two 
shifted constants only. 
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2.3 Multiscale Representation 



Accordingly to our approach which allows us to find exact solutions as for 
Vlasov-like systems as for rms-like systems we need not to fix particular case of 
distribution function F(x, y, x' , y', s). Our consideration is based on the follow- 
ing multiscale TV-mode anzatz: 

N 5 

F N (x,y,x',y',s)= ]T a iu ... ii5 <^}A ik (x,y,x',y',s), (30) 

»i,...,»5 = l k=l 
N 3 

i/j N (x,y,s) = & iij2,i 3 0%(^.2/>s). (31) 

31,32,33 = 1 fc=l 

These formulae provide multircsolution representation for variational solutions 
of our system. Each high- localized mode/harmonics Aj(s) corresponds to level 
j of resolution from the whole underlying infinite scale of spaces: 

...V-2 C V-i C V C Vi C V 2 C (32) 

where the closed subspace Vj(j £ Z) corresponds to level j of resolution, or 
to scale j. The construction of tensor algebra based on the multiscale bases 
is considered in part I [7]. We will consider rate equations as the following 
operator equation. Let L, P, Q be an arbitrary nonlinear (rational in dynamical 
variables) first-order matrix differential operators with matrix dimension (d = 4 
in our case) corresponding to the system of equations, which act on some set of 
functions 

* = *(«) = (V(s),...,* d (s)J, sencR (33) 

from 

L 2 (fl) : Q(R, «)*(«) = P(R, s)*(s) (34) 

or 

=L(fl,s)tf(s) = (35) 

where 

R = R(s,d/ds,^>). (36) 

Let us consider now the N mode approximation for solution as the following 
expansion in some high-localized wavelet-like basis: 

N 

* N (s) = J2aUr(s). (37) 

r=l 

We will determine the coefficients of expansion from the following variational 
condition: 

L% = J (M N )<j> k (s)ds - (38) 

We have exactly dN algebraic equations for dN unknowns a r . So, variational 
approach reduced the initial problem to the problem of solution of functional 
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equations at the first stage and some algebraic problems at the second stage. 
As a result we have the following reduced algebraic system of equations (RSAE) 
on the set of unknown coefficients af of the expansion: 

HiQi^a?,^) = MiP^afrfj), (39) 

where operators H and M are algebraization of RHS and LHS of initial problem. 
Qij (Pij) are the coefficients of LHS (RHS) of the initial system of differential 
equations and as consequence are coefficients of RSAE (39). 

I = (ii,...,i q+2 ), J = (ji, j p +i) (40) 

are multiindices, by which are labelled otj and /?/, the other coefficients of RSAE: 

Pj = {Pn...j P+1 }= f LI ( 41 ) 

l<jk<p+l 

where p is the degree of polynomial operator 

ai = {a^...a iq+2 } = J ^...^...cf)^, (42) 

il,...,iq+2 

where q is the degree of polynomial operator Q, 

U = (l,...,q + 2), <4=d^ s /ds. (43) 

We may extend our approach to the case when we have additional constraints 
on the set of our dynamical variables 

* = {i?,e} (44) 

and additional averaged terms also. In this case by using the method of La- 
grangian multipliers we again may apply the same approach but for the extended 
set of variables. As a result we receive the expanded system of algebraic equa- 
tions analogous to our system. Then, after reduction we again can extract from 
its solution the coefficients of the expansion. It should be noted that if we con- 
sider only truncated expansion with N terms then we have the system of N x d 
algebraic equations with the degree 

I = max{p, q} (45) 

and the degree of this algebraic system coincides with the degree of the initial 
system. So, after all we have the solution of the initial nonlinear (rational) 
problem in the form 

N N 

R N ( S ) = R(0) + J2aUk(s), e N (s) = e(0) + (46) 
fe=i fe=i 
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where coefficients 

a%, b» (47) 

are the roots of the corresponding reduced algebraical (polynomial) problem 
RSAE. Consequently, we have a parametrization of the solution of the initial 
problem by solution of reduced algebraic problem. The problem of computations 
of coefficients 

«/, Pj (48) 

of reduced algebraical system may be explicitly solved in wavelet approach. 
The obtained solutions are given in the form, where </>fc(s) are proper wavelet 
bases functions (e.g., periodic or boundary). It should be noted that such rep- 
resentations give the best possible localization properties in the corresponding 
(phase)space/time coordinates. In contrast with different approaches these for- 
mulae do not use perturbation technique or linearization procedures and repre- 
sent dynamics via generalized nonlinear localized eigenmodes expansion. Our 
N mode construction gives the following general multiscale representation: 

R{s) = <°™( S ) + ^ff(^ s ), uj^T, (49) 

i>N 

e(s) = e#°»+£y( Wj -s), u^2\ 

j>N 

where R l (s), e J (s) are represented by some family of (nonlinear) eigenmodes 
and gives the full multiresolution/multiscale representation in the high-localized 
wavelet bases. As a result we can construct various unstable (Fig. 10, Part I) 
or stable patterns (Fig. 1, 2 , Part II and Fig. 11, Part I) from high-localized 
(coherent) fundamental modes (Fig. 1 or 8 or 9 from Part I) in complicated 
stochastic systems with complex collective behaviour. Definitely, partition(s) 
as generic dynamical variables cannot be postulated but need to be computed 
as solutions of proper stochastic dynamical evolution model. Only after that, 
it is possible to calculate other dynamical quantities and physically interesting 
averages. 



3 CONCLUSIONS: TOWARDS ENERGY 
CONFINEMENT 

Analysis and modeling considered in these two Parts describes, in principle, 
a scenario for the generation of controllable localized (meta) stable fusion-like 
state (Fig. 2 and Fig. 6). Definitely, chaotic-like unstable partitions/states 
dominate during non-equilibrium evolution. It means that (possible) localized 
(meta) stable partitions have measure equal to zero a.c. on the full space of 
hierarchy of partitions defined on a domain of the definition in the whole phase 
space. Nevertheless, our Generalized Dispersion Relations, like (15) or (39), 
give some chance to build the controllable localized state (Fig. 6) starting from 
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initial chaotic- like partition (Fig. 3) via process of controllable self-organization. 
Figures 4 and 5 demonstrate two subsequent steps towards creation important 
fusion or confinement state, Fig. 6, which can be characterized by the presence 
of a few important physical modes only in contrast with the opposite, chaotic- 
like state, Fig. 3, described by infinite number of important modes. Of course, 
such confinement states, characterized by zero measure and minimum entropy, 
can be only metastable. But these long-living fluctuations can and must be very 
important from the practical point of view, because the averaged time of exis- 
tence of such states may be even more than needed for practical realization, e.g., 
in controllable fusion processes. Further details will be considered elsewhere. 
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Figure 1: Localized two-dimensional partition. 




o o 

Figure 2: Metastable pattern. 
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Figure 3: Chaotic pattern. 




o o 

Figure 4: Self-organization, step 1: towards confinement. 
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Figure 5: Self-organization, step 2: towards confinement. 




o o 

Figure 6: Localized pattern(waveleton): energy confinement state. 
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